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SUPPORT THEOREMS FOR THE LIGHT RAY TRANSFORM 
ON ANALYTIC LORENTZIAN MANIFOLDS 

PLAMEN STEFANOV 


Abstract. We study the weighted ray transform of integrating functions on a Lorentzian manifold 
over lightlike geodesics. We prove support theorems if the manifold and the weight are analytic. 


1. Introduction 

Let g be a Lorentzian metric with signature on the manifold M of dimension 

1 + n, n > 2. Light-like geodesics 7 (s) (also called null geodesics) are the solutions of the geodesic 
equation = 0 for which g( 7 , 7 ) = 0. There is no canonical unit speed parameterization as in 
the Riemannian case as discussed below. For some fixed choice of it, we define the weighted light 
ray transform L K f of a function (or a distribution) / on M by 

( 1 - 1 ) L K f(i) = J «( 7 (s), 7 (s))/( 7 (s))ds, 

where 7 runs over all null geodesics. Here k is a weight function, positively homogeneous in its 
second variable of degree zero, which makes it parameterization independent. When k = 1, we use 
the notation L. Conditions for supp / and the interval of definition of the geodesics will be specified 
below but in all cases, the integration is carried over a compact interval. This transform appears 
in the study of hyperbolic equations when we want to recover a potential term from boundary or 
scattering information, see, e.g., EH 123 121 EH m EH M for time dependent coefficients, and also 
|2) [22] for time-independent ones. It belongs to the class of the restricted X-ray transforms since 
the complex of curves is restricted to the lower dimension manifold < 7 ( 7 , 7 ) = 0. 

Our goal is to study the invertibility of L K , including its microlocal invertibility. While the 
methods we develop could be used to study stable recovery of the ( C °°) spacelike wave front set 
of /, we concentrate our attention here on support theorems for analytic metrics and weights. In 
[29] . the author showed that if g is the Minkowski metric, and if f(t,x) is supported in a cylinder 
R x B(0,R) and has tempered growth in the time variable, then Lf determines / uniquely, see 
also [25] . The proof was based on the fact that Lf recovers the Fourier transform / of / (w.r.t. all 
variables) in the spacelike cone |t| < |£| in a direct (and stable) way and since /(r, £) is analytic in 
the £ variable (with values distributions in the t variable), then one can fill in the missing cone by 
analytic continuation in the £ variable. It is easy to see that there is no stable way to recover / in 
the timelike cone |r| > |£| (true also in the most general Lorentzian case, see next paragraph) thus 
L has a high degree of instability, see also [T]. From a physical point of view, this could be expected: 
we can recover all “signals” moving slower than light, and we should not expect to recover those 
moving faster than light; and the latter should not exist anyway expect for possible group velocity 
faster than light. 
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When the metric is not flat, it is fairly obvious that L K f cannot “see” the wave front WF(/) 
in the timelike cone, this just follows from the inspection of the wave front of the Schwartz kernel 
of L k , see also Theorem m for the Minkowski case. Recovery of WF(/) in the spacelike cone 
is far less obvious and certainly requires some geometric assumptions like no conjugate points or 
existence of a foliation of strictly convex surfaces, as we explain below. One possible approach is to 
analyze the normal operator L' K L K as in mmm- That operator is in the I P)l class of 'I'DOs with 
singular kernels, which are Fourier Integral Operators (FIOs), in fact, see [lOj and the references 
there. The analysis of L' K L K in the Minkowski case for n = 2 is presented in ME] as an 
example illustrating a much more general theory. Applying the I p ’ 1 calculus to get more refined 
microlocal results however requires the cone condition which cannot be expected to hold on general 
Riemannian manifolds due to the lack of symmetry, as pointed out in m- An alternative approach 
to recover the C°° spacelike singularities can be found in m- 

Our main result is support theorems and injectivity of L K for analytic metrics and weights (on 
analytic manifolds M). It can be viewed as an extension of the classical Helgason support theorem 
for Radon transforms in the Euclidean space pH]. We use analytic microlocal arguments. Such 
techniques go back to MU- In [5], the authors prove support theorems for Radon transforms 
(with flat geometry) and analytic weights. In [6], they study “admissible line complexes” in R 1+2 
with analytic weights, and type III there includes a weighted version of L K in the Minkowski case. 
Their arguments however are based on the calculus of the analytic FIOs as an analytic version of 
the C°° analysis in [8]. Such a generalization does not exist to the best of the author’s knowledge. 
Even the analytic TDO calculus is quite delicate already, see, e.g., [35], and an analytic version of 
the FIO calculus, including the I p,l one, would pose even more challenges. Support theorems for 
the geodesic transforms on simple analytic manifolds have been proved with analytic microlocal 
techniques in [E]CE| and related results; even for tensor fields in |7][32][33]- A breakthrough was 
made by Uhlmann and Vasy in [36]; who proved a support theorem in the Riemannian case near 
a strictly convex point of a surface in dimensions n > 3 without the analyticity condition. The 
X-ray transform is assumed to be zero on all geodesics close to tangent ones to the surface at that 
point, and / is a priori supported on the concave side. Their arguments are based on application 
of the scattering calculus [20] and the n > 3 assumption is needed to guarantee ellipticity in a 
neighborhood of the point. 

The approach we propose is simpler and avoids all the difficulties related to the singularities of 
the symbol of L' K L K : we form smooth timelike surfaces foliated by lightlike geodesics over which 
one can compute a weighted Radon transform R by just applying Fubini’s theorem. This reduces 
the problem to a microlocal inversion of that (non-restricted) Radon transform known on an open 
set of surfaces, which in the smooth case is doable with classical microlocal techniques going back 
to Guillemin nans]. Analytic singularities can be resolved by the local Radon transform, as well 
m • On the other hand, this approach does not allow us to analyze the lightlike singularities, where 
some form of the I p ’ 1 calculus would still be needed. In the proof of Theorem 14. II those surfaces do 
not appear explicitly but they can be thought of as the level surfaces of the phase function q i. One 
would expect to be able to do the analytic microlocal inversion by treating R'R as an analytic TDO 
but it is not clear how to do that to obtain purely local results due to the delicate nature of cut-offs 
allowed in that calculus. Instead, we use the analytic stationary phase approach by Sjostrand [28] 
already used by the author and Uhlmann in [33], see also Humus]. 

As a simple example illustrating the reduction of the restricted ray transform L K to a classical 
Radon transform R , consider the Minkowski case. Light geodesics are given then by the lines parallel 
to (1 ,8), with |@| = 1. Every timelike plane (with a normal v = (^, v x ) such that \vt\ < |i/ x |) can 
be represented easily as a foliation of light rays. If L K f € C°° (or analytic), then so is Rf on the 
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open manifold of those planes. Then we have to invert microlocally the classical Radon transform, 
which is well known. This argument still works if we introduce a weight in L K and/or know L K f 
localized to an open set of light rays only, see Theorem 13.21 

Finally, we notice that some global conditions on the geodesic flow are clearly needed for mi¬ 
crolocal inversion, even in the spacelike cone. If g = —At 2 + hij(x) dxMx- 5 , where h is a Riemannian 
metric on a bounded domain, then L K , restricted to f-independent function reduces to the geodesic 
X-ray transform X. It has been shown recently mm that when n = 2, Xf recovers WF(/) in 
a stable way if and only if there are no conjugate points. When n > 3, the no conjugate points 
condition is sufficient [32, 33] and there are examples of metrics of product type for which it is 
necessary, by the 2D results in pH [21] . On the other hand, the support theorem in [36] provides 
global uniqueness and stability under another type of condition: existence of a foliation by strictly 
convex surfaces (conjugate points may exist). This implies stable invertibility when n> 3 without 
analyticity assumptions. We assume the foliation condition but in contrast to [36] . here n = 2 
is allowed since for our purposes, full ellipticity (in all directions at a point) is not needed; only 
ellipticity at directions conormal to the foliation suffices. Also, full ellipticity does not hold in the 
Lorentzian case since the timelike singularities are invisible. On the other hand, we require g and 
k to be analytic. One would expect support theorems under the no-conjugate points assumption 
as well but that remains an open question. 

Acknowledgments. The author thanks Manuel Gutierrez for helpful discussions about Lorentzian 
geometry and to Gunther Uhlmann for numerous discussions about Integral Geometry and Inverse 
Problems. 


2. Main results 

2.1. Support theorems for the Minkowski spacetime. Let g = — df 2 + (dx 1 ) 2 -|-b(dx n ) 2 be 

the Minkowski metric in R 1+Tl . Future pointing lightlike geodesics are given by s H > (t + s, x + s6) 
with |0| = 1. They can be reparameterized by shifting and rescaling s. Note that the notion of 
“unit” speed is not invariantly defined under Lorentzian transformations but in a fixed coordinate 
system, the scaling parameter 1 (i.e., At/ds = 1) is a convenient choice. Set 

(2.1) Lf(x,0) = J f(s,x + sO)As, xeR n ,6eS n ~ 1 . 

This definition is based on parameterization of the lightlike geodesics (lines) by their point of 
intersection with t = 0 and direction (1,0). We will use the notation 

(2.2) 4,6 i(s) = (s,x + s9). 

The parameterization (x, 6) defines a natural topology and a manifold structure of the set of the 
future pointing lightlike geodesics. 

Given a weight k € C'°°(R x R” x S'” -1 ), we can define the weighted version L K of L by 

L K f(x,0) = j k(s,x + s0,9)f(s,x T sO) ds, x € R”, 0 € S n-1 . 

In the terminology of relativity theory, vectors v = (v°,v') satisfying |uo| < |u'| (i.e., g(v,v) > 0) 
are called spacelike. The simplest example are vectors (0,u'), v' ^ 0. Vectors with |uq| > |u'| (i.e., 
g(v,v) < 0) are timelike; an example is (1,0) which points along the time axis. Lightlike vectors 
are those for which we have equality: g(v, v ) = 0. For covectors, the definition is the same but 
we replace g by g , which is consistent with the operation of raising and lowering the indices. 
Surfaces with timelike normals (which are covectors) are spacelike, and vice versa. 
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Definition 2.1. Let K be a subset of the Minkowski spacetime. We say that I\ expands with speed 
less than one if 

(2.3) K C \x\ < c\t\ + R} for some 0 < c < 1, R > 0. 

Condition (12.311 is easily seen to be invariant under Lorentz transformations. Also, it does not 
require supp / to be compact. 

The terminology we used is a bit ambiguous. What we actually mean is that the cross-section 
of K with any plane t = const, is a bounded set contained in a ball expanding with a speed less 
than one. If dK is smooth, we do not really require it to be timelike. 

In the Minkowski spacetime, we have the following support theorem. 

Theorem 2.1. Let f E V'(Yl l+n ) be so that supp / expands with a speed less than one. Let l X0 ,e 0 
be a fixed lightlike line in the Minkowski spacetime and let U 3 (xo,0o) be on open and connected 
subset of R n x S n ~ 1 . Let n(t,x,0) be analytic and non-vanishing for (t,x) near supp / so that 
(x — tO , 0) E U. 

If L K f(x, 0) = 0 inU and if I X0 ,6 0 does not intersect supp f, then none of the lines I Xi q, (,r, 0) E U, 
does. 

2.2. Support theorems on Lorentzian manifolds. Let ( M,g ) be a Lorentzian manifold. Light¬ 
like (null) geodesics are defined as the geodesics 7(5) for which <7(7,7) = 0. They exist at least 
locally by the ODE theory. There is no canonical parameterization since for any linear transfor¬ 
mation of the s variable cr(s) = as + b, a 7 ^ 0, 7 o a is still a null geodesic. Moreover, a and b may 
change from geodesic to geodesic. Let S be a spacelike surface near a fixed lightlike geodesic 70 (s), 
intersecting S for s = 0. Then we can parameterize the lightlike geodesics in some neighborhood of 
7o(0) E S close to 70 with directions close to 7(0) with initial points 2; on S' and initial lightlike di¬ 
rections v at x pointing in the direction of 70. A choice of the scaling of the parameter s along each 
7(s) can be fixed by requiring 7(0) E S and requiring the normal component of 7 on S to be a given 
negative function, for example —1. If that function is smooth/analytic when S is smooth/analytic, 
we call the parameterization smooth/analytic. This property does not depend on the choice of S 
and also defines a topology and a smooth/analytic structure of the lightlike geodesics defined on a 
fixed interval. We could use a timelike surface as initial points instead. 

If C C M is closed, we call the null geodesic y(s) non-trapping in C, if 7 -1 (C) is contained in some 
open finite interval call it / for a moment. For any local parameterization of null geodesic as above, 
the maximally extended null geodesic with initial points and directions close enough to 7 would 
leave C for s near the ends of /. Some of them may return to C for s 0 / (even though this cannot 
happen to 7 but we restrict them to / only. Then we consider those geodesics a neighborhood of 
7, identified with the neighborhood of the initial points and directions in that parameterization. 
That definition of local neighborhood is independent of the chosen parameterization and defines a 
topology near 7 (restricted to I). For any such choice of the parameterization, we then define L K f 
locally by m for any / E with s restricted to I. A different analytic parameterization would 
change L K f (in a trivial way) but it will not change its property to be smooth or analytic, or zero. 

We do not want to assume that / is compactly supported but we always assume that we integrate 
over a set of light geodesics non-trapping in supp/. Then locally, we may cut, a smooth / in a 
smooth way to make it compactly supported without changing L K f near that geodesic. This 
reduces the local analysis to compactly supported functions. An example is a function supported 
in a cylinder \x\ < R in the Minkowski case; or more general / with supp / expanding with speed 
less than one, see (12.31) . This allows L K f to be well defined for smooth / over open sets of 11011- 
trapped light geodesics and then by duality for distributions /. Indeed, for every distribution / 
in M we can set locally L n f = L K \f , with a suitable x € Cq°, and the latter makes sense by 
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duality. In other words, near a fixed null geodesic, non-trapping in supp /, it is enough to study 
L k restricted to compactly supported distributions /. Based on that, to simplify the formulation 
of the next theorem, we assume that supp / is compact. 

Theorem 2.2. Let (M, g ) be an analytic Lorentizan manifold and let k be an analytic non-vanishing 
weiqht. Let F : M —>• [0, ll be a smooth function. Assume that f £ £'(M) and 

(i) F _ 1 ( 0 ) fl supp / = 0 , 

(ii) d F / 0 on supp /, 

(in) F^ 1 (a) fisupp/ is strictly lightlike-convex for all a £ [ 0 , 1 ]. 

Then if L K f( 7 ) = 0 in a neighborhood of all null-geodesics with the property that each one is tangent 
to some of the surfaces F~ 1 (a), a £ [0,1], then f = 0 on T -1 !/), 1). 

We refer to Definition 14.11 for the notion of lightlike convexity. Examples of strictly lightlike- 
convex surfaces in the Minkowski spacetime, which cannot be spacelike at any point, include the 
cylinder |x| = R, R > 0; more generally, the smooth part of the double cone \x\ = c\t\, with 
0 < c < 1 fixed; or the hyperboloid \x\ 2 = eft 2 + C with C > 0 and such a c. They are all timelike. 
We also note that we can actually require L K f = 0 on a suitable submanifold of lightlike geodesics 
of dimension 1 + n only, as it follows form the proof. Moreover, it is only enough to assume that 
L K f is analytic there, even microlocally so, see Remark 13.11 

To demonstrate a typical application of Theorem 12.21 we will point out how one can show that 
if in the Minkowski space time / satisfies (12.31) . and Lf = 0, then / = 0, which, of course follows 
from Theorem 12.II as well. We choose F(a) = a + c 2 (t — to) 2 — \x — xo\ 2 with 0 < c < c < 1, <7 > 0. 
The constant a can be rescaled to be fit in [0,1]. Choosing various xo and to we can prove / = 0. 
One can perturb the metric a little bit assuming that supp / is supported in a fixed compact set, 
to get examples for non-flat metrics. 

3. Analysis in the Minkowski case 

3.1. Fourier Transform analysis. Let M = R” and g be Minkowski. By the Fourier Slice 
Theorem, knowing the X-ray transform for some direction oj recovers uniquely f on oj 2 -. More 
precisely, the Fourier Slice Theorem in our case can be written as 

/W*=o = f(-0 ■£,£) = [ e~ ix <Lf(x, 6) dx, V9 £ 5”' 1 . 

Here and below, we denote by ( = (r, £) the dual variables to z = (t,x). The proof is easy, see 
(13.111) . The union of all (1,0) _L for all unit 6 is {|r| < |£|} = U L/, that is easy to see. This 
correlates well with the theorems below. In particular, we see that knowing /(£) for a distribution 
/ with a well defined Fourier transform, recovers / in the spacelike cone uniquely and in a stable 
way. Under the assumption that supp / is contained in the cylinder \x\ < R for some R (and 
temperate w.r.t. t), one can use the analyticity of the partial Fourier transform of / w.r.t. x to 
extend / analytically to the timelike cone, as well. This is how it has been shown in [29] that L is 
injective on such /. 

3.2. The normal operator X'X. We formulate here a theorem about the Schwartz kernel of the 
normal operator N = L'L. We will skip the proof because we will not use the theorem for our main 
results. One way to obtain it is to think of L as a weighted version of the X-ray transform X with 
a distributional weight 5(t 2 — |£| 2 ) and use the results about the weighted X-ray transform, see, 
e.g., m and allow a singular weight there. Details will appear in [30], see also [19] . 
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Theorem 3.1. 

(a) 


L'Lf = Af * f, Af(t,x) = 


5(t — |x|) + 5(t + |x| 

\rr\n-l 


(b) 


(c) 


n —3 

( 1£|2 _ t 2 \ 2 

L'Lf = C n F- 1 " q \„ n J + Ff, V/ e S(R 1+ri ), C n :=2n\S' 


n— 2 1 


Ifl 


3 —n 

h(0 + )f = C- 1 \D x \ n ~ 2 D + 2 X'X /, 


where h is the Heaviside function, and □ = d? — A z and J- is the Fourier transform. 


Above, we used the notation = max(s, 0) with the convention that is the Heaviside 
function. 

In particular, when n = 3, we get a(X'X) = Cs\^ 1 h (|£| 2 — |£o| 2 )- Then 

/»(□+)/ = C^\D z \X'Xf. 

As we can expect, there is a conormal singularity of the symbol even away from £ = 0 living on 
the characteristic cone, and X'X is elliptic outside it, and only there. The theorem shows that 
“singularities traveling slower than light” can be recovered. The ones traveling faster cannot. 


3.3. Recovery of spacelike C°° singularities in the Minkowski case. Our first theorem says 
that knowing Lf near a lightlike geodesic £q allows us to recover all spacelike singularities conormal 
to £q. We denote the conormal bundle of £q by N*£q. Recall that the conormal singularities to 
£q contain lightlike ones, as well. This result follows from the analysis in nmm and the reason 
we present it here is to illustrate the main idea on a simpler problem where we can do explicit 
computations. 

Theorem 3.2. Let f be a distribution so that l(s) ^ supp / for |s| > 1/C with some C for all 
lightlike lines £ near £$. Let L K f(£) € C°° for £ in some neighborhood T of £q. Then WF(/) n N*£q 
contains no spacelike covectors. 

Proof. We construct planes foliated by lightlike geodesics with a fixed direction (1, 6) € R x S n ~ l . 
Any such plane intersects the t = 0 plane in a (n — l)-dimensional plane in the x space. Let the 
latter be ir P)U} = {x • u = p}, u € 5 n_1 , p G R. Then the plane that we denote by n is the flow 
out of the null geodesics with celestial direction 6 originating from 7r PA ,, i.e., 

(3.1) ir PtUJ ,e = {(t, x + td)\ x € 7T PtUJ }. 

The same plane can be also described by the equation (x — tO) ■ u = p, therefore, 

(3.2) iT PjW} g = {(t,x); (t,x) ■ (-0 ■ u,u) = p}. 

The dot product here is in the Euclidean sense, and can be also thought of as a pairing of a vector 
and a covector, in invariant terms. In particular, we see that the set of such planes coincides with 
the timelike ones and the lightlike ones as a borderline case. 

Let C° ^ 0 be spacelike and conormal to £$ at a point that we can always assume that to 
be the origin. Applying a Lorentz transformation, we can always assume that = e n_1 := 
(0,..., 0,1, 0) € R 1+n and £o = £o,e„ = (s, 0,... 0, s). Here and below, we use the notations 
and e k to denote vectors/covectors with all entries zero instead of the k -th one. Take the plane 
7 T0 = {(i, x)', x n ~ L = 0}, conormal to £°. This is the plane constructed above with u = £° = e 
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and 8 = e n (we could have chosen any other 8 _L u> but we chose this one because it is related to 


^o)- 

It is more convenient to extend the parameters (p, uj) by homogeneity. We allow uj to be non-unit 
and denote it by £. Then the planes 7r ' PjU g are given by 


(3.3) z-C = P, ( = (-8-Z,0- 

We will choose a suitable analytic family of (£, 0) near £ = e n_1 , 9 = e n parameterized by an 
n + 1-dimensional parameter so that the map from that parameter to the (co-)normal £ is a local 
diffeomorphism. We keep £ unrestrained and let 9 depend on an ID parameter, q. Then 


(3.4) dC/d^k = (~8 k ,e k ), k = l,...,n, and dC/dq\ q=0 = {-89/dq ■ uj\ q=0 , 0). 


This system of vectors is linearly independent, if and only if d9/dq ■ uj | g= o / 0. Therefore, the 
variation 89/dq should be chosen not parallel to 7T p ^. This leaves essentially a variation in the 
direction of £. Based on that, we set 


(3.5) 9(q) = (cos q)e n + (sing)e n _i. 

Then 


(3-6) £(<?,£) = (-%) ■£,£) 

and 


(3.7) 7Tp^,0(g) = {{t, x); (t, x ) • C (q, C) =p}i P e R; £ g R n \ 0. 

The fact that (g, £) —> £ is a local diffeomorphism is also easy to verify directly. Solving (13.61) for 
(q, £) yields £j = Ci, i ~ 1, ■ ■ ■, n, and 

(3.8) Cn cos q + £ n _i sin q = —£ 0 

and the latter is uniquely solvable for q near q = 0 for £ near e n ~ 1 ; let q = g(£) be the solution. 

We can write the defining equation also as (t,x) • v = p := p/\(\ with v = £(g, £)/|<£(<3', £)|- Then 
it is easy to show that with £ restricted back to unit sphere, the map R x S n ~ l 3 (g,£) v £ S n 
is a local analytic diffeomorphism near p = 0. In other words, (q. £) with £ unit, parameterizes the 
normal v to (13.71) in a locally diffeomorphic way. 

By the support assumption of the theorem, there exists a neighborhood U of (0, e n ) and A > 0, 
so that all lightlike geodesic issued from U leave supp/ for |s| > A. Take a smooth function x( x ,9) 
supported in U equal to 1 near (0,e n ). Since L K f G C°°, we have X L Kf\xen p ^ gM ,O=0(q) € C°°, as 
well. Integrate [xL K f](x,6(q)) with respect to x on the plane n p ^ to get by Fubini’s theorem: 

(3-9) Rf(n p ^ e{q) ) := f X«/djU Pl g,o(,) G C°° 

£,«(«) 

for some measure analytically depending on (p, q, £), i.e. an analytic and positive multiple of the 
Euclidean measure on each plane. Above, the integral is taken over a compact set; moreover, we 
can cut / to a compactly supported distribution away from where we integrate without affecting 
the integral. Therefore, / is in the microlocal kernel of the weighted Radon transform R with a 
weight not vanishing at (t,x) = 0 on the plane to- This allows us to apply R' to get an elliptic 
\kDO of order —2, see, e.g., [6]. Therefore, / is microlocally smooth at (0,£°) as claimed. □ 

Remark 3.1. Note that we only needed to know that L K f(x,9) vanishes (being microlocally 
smooth in a certain cone would suffice) for 9 = 9(q) only with |g| <C 1; i.e., we require knowledge 
of a restricted version of the already restricted L. This is similar to the known fact that in the 
Euclidean space we can invert the X-ray transform by “slicing” R n into 2D planes. We could have 


P. STEFANOV 


proven our results in the Minkowski spacetime in 1 + 2 dimensions only and then extended it to 
any dimension 1 + n > 1 + 2. The same remarks applies to the analytic case below but we need 
to know that L K f is microlocally analytic (instead of just smooth) in some conic set. Even in the 
Lorentzian case, we still need to know L K restricted to a certain an (l + n)-dimensional submanifold 
of geodesics. 


3.4. Recovery of analytic spacelike singularities in the Minkowski case. We show first 
that we can recover all spacelike analytic singularities of / conormal to the lightlike lines along we 
integrate. For a definition of the analytic wave front set WFa(/), we refer to [28] and [35] . 

Lemma 3.1. Let f G P / (R 1+ri ) and let t X0 ,e 0 be a fixed lightlike line so that £ x ,e{s) fL supp / for 
|s| > 1/C with some C for all (x,6) near (xo,0o)- Let k(s,x + s9) be analytic and non-vanishing for 
those s, x, 0. If L K f(x, 6) = 0 near (xq, 9q), then N*l X0t g 0 n WFa(/) contains no spacelike covectors. 


Proof. One would expect the proof to be a complete analog of that of Theorem 13.21 but that proof 
involves smooth cutoffs along the planes we integrate over. We cannot do this in our case because 
that would destroy the analyticity of the weight. On the other hand, we need the localization 
because we know that L K f = 0 near I XQ ,0 o only. 

We use the local coordinates in the proof of Theorem 13.21 where C° = e n ~ 1 , 9 q = e n , zq = 0. 

Let xn € C'^°(R” ) be with support in B(0, e), e > 0, with xn = 1 near xq = 0 so that 

(3.10) \9 x Xn\ < (CN)\ a \ for |a| < N, 

see [35]. Then for 0 < e <C 1, A > 0, and 0 close to 9q = e n , 


0 = j e lXx< ( X NL K f)(x,0)dx 


e lXx ^X n(x)k(s, x + s9, 9)f(s, x + sQ) ds dx. 


If (1, 6) ■ f = 0 with f = (r, f), then x ■ £ = (t, x + t6) ■ f. Make the change of variables x + s9 <—> x 
in the integral above to get 


(3.11) 


0 = 


e lXx< (xNL K f)(x,9)dx 


= fj e iA(t ' lK Xiv( 


x — t6)n(t , x , 9)f(t, x) df dx, 


if (1,0) -c = o. 


For ^ G R n , choose 9 = 9{q) as in ( 13 . 51 ) . and set £ = (— 9(q) -f, as in ( 13 . 31 ) . Then the orthogonality 
condition in ( 13 . 111 ) is satisfied. To connect this with the analysis in section [3731 notice that we can 
get the same result by taking the Fourier transform IF p ^,\ in ( 13 . 91 ) . with y = xn- Choose now 
q = q(Q as in ( 13 . 81 ) . The orthogonality condition still holds and we have 


(3.12) 


J e lA2 '^a 7 v(+ ()f(z) dz = 0 near £ = e n ^ 


where ajv = Xn(x — t9(q))n(t, x, 6(q)), with q = q {£), is analytic and elliptic near (z, Q = ( 0 , e n ~ l ) 
(but not analytic away from some neighborhood of it) and satisfies pseudo-analytic estimates of 
the type ( 13 . 101 ) . 

We will apply the complex stationary phase method of Sjostrand [28] similarly to the way it was 
applied in m to the partial data Calderon problem and in mm to integral geometry ones. 

Fix 0 < 4 <C £, see ( 13 . 101 ) . With some w € R 1+n , r? G R 1+n close to w = 0 , n = e n_1 , multiply 
the l.h.s. of (l+TTl) by 


( 3 . 13 ) 


e iA(i(C-+ 2 /2-«K) 
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and integrate w.r.t. ( in the ball \C ~ v\ < <5 to get 


(3.14) 



e iX ^ w ’ z ^ aN ( Z ,C)f( Z )dzd( = 0, 


where 

<t> = (z — w) ■ ( + i(( — r\) 2 /2 . 


We split the z integral (13.1411 into two parts: over {z; \z — w\ < 5/2} and then over the complement 
of that set. Since |4>^| has a (5-dependent) positive lower bound for \z — w\ > 5/2, we can integrate 
in the outer integral in (13.1411 by parts w.r.t. (, see, e.g., [7] [33] using (13.1011 and the fact that on 
the boundary |( — rf\ = 5, the factor e lA$ is exponentially small with A. We then get 


(3.15) 



l«5/2, |C-??|<<5 


e iX ^ w ’ z ’^ a { Z ,C)f{z) dzdC 


< C{CN/X) N + CNe~ x/c 


where a equals a with the \N factor missing, i.e., a = n(t,x,9(q(Q), which is independent of N 
because on the support of the integrand, that factor is equal to 1 for 5 <C e, see (|3.12p . Choose 
now N so that N < \/{Ce) < N + 1 to get an exponential error on the right. 

The phase <f>, as a function of (, has a unique critical point ( c = rj + i(z — w) and |C c — r}\ < 5/2 
on the support of the integrand in (13.151) . Set 


(3.16) 


il>{w,z,rj) = $| C =Cc- 


Therefore, 

ip = f] ■ (z — w) + i\z — w\ 2 — —\z — w\ 2 = rj ■ (z — w) + -\z — w\ 2 . 

This is the type of phase functions that are used to test for analytic microlocal regularity. We apply 
now [28, Theorem 2.8] and the remark after it to the ("-integral in ([3.151) to get 


(3.17) 


'\z—w\<5/2 


e iX ^ w ’ z ^ b (w,z, r i,X)f{z) dz 


< Ce~ x ' c . 


for (z,r]) close to (0,e n_1 ), with some classical elliptic analytic symbol b of order 0 in the sense of 
[28] near (w, z, rj) = (0,0, e n_1 ). In particular, the principal part of 6(0,0, e n , A) is f3n( 0, 0, e n ) with 
/3 an elliptic factor depending on the phase, see [28], Theorem 2.8]. This implies (0, e"^ 1 ) 0 WFa(/), 
see [28] ■ □ 


3.5. Proofs of the support theorems in the Minkowski spacetime. The next proposition 
is a unique continuation result across a timelike surface in the Minkowski case which implies The¬ 
orem m 

Proposition 3.1. Let f E P / (R 1+n ) and let £ X0) e 0 6e a fixed lightlike line in the Minkowski space- 
time so that £ x ,e(s) 0 supp / for |s| > 1/C with some C for all {x, 9) near (xo, 9 q)- Let k(s, x + sQ) 
be analytic for those s, x, 9. 

If L k J(x , 9) = 0 near (®o, ^o) if f = 0 on one side of S near zo, then f = 0 near zo. 

Proof. Assume that zo € supp/. Then (zo, tK-o)) € WFa(/) by the Sato-Kawai-Kashiwara 
Theorem, see m and [28], where v(zq) is one of the two unit co-normals to S at zo- That covector 
is spacelike by the assumption about S, and is conormal to £ Z0 p 0 (ff). This contradicts Lemma [Til 
which completes the proof of the proposition. □ 
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Proof of Theorem \2.1[ Fix {x\,9\) E U. Let [0,1] 3 p —» (x p ,9 p ) be a continuous family in U 
connecting (xo,#o) with ( 21 , 6 b). We can always assume that U is bounded, hence U is compact. 
Let U be the set of points lying on £ Xj g, ( 2 , 9) E U. 

Choose c E (c, 1), where c is the constant in (|2.3j) . Denote by £ Xt g(s) = (s,x + sc9), (x,9) E 
R n xS n 1 , the timelike geodesics with speed c. By (12.31) . there exists A > 0 so that £ Xj g(s) & supp / 
for all (x,6) E U and |s| > A, and so that the same holds for £ X) q{s ) uniformly w.r.t. c as long as 
c + n < c < 1 with some fixed /iE (0,1 — c). 

Let e > 0 be such that the cylinder C X0t g 0 := Ui y - Xo \< £ £y,e is disjoint from supp / but sill lies in 
U. By the arguments above, we can assume that supp / is compact since t is bounded on it along 
the light lines under consideration. Therefore, for the cylinder C x g := ^\ y - x \<E./ 2 ^y,e we have 

(3.18) C x>e n {\t\ <A}c C X} e n {\t\ < A} 

for every ( x , 9) as long as c is close enough to c but still smaller than it. We require one more 
property for c which actually refines (|3.18l) : 

V(x, 9) E U,\/z E C X) e D {11| < A}, and every unit 9\ with \9\ — 9\ < yjl — c 2 , 
the lightlike line through z in the direction of (1, 9 \) stays in C x> q D [\t] < A}. 

This property can be guaranteed for 1 — c <C 1 by continuity and compactness. We fix such a c. 

Assume that the family {C x q ; p E [0,1]} has a common point with supp/. Let po be the least 
p E [0,1] (which exists by compactness and continuity arguments) for which C Xp j) p H supp / ^ 0. 
Then / = 0 in the interior of C Xp .o p and there is a point on its boundary which is also in 
supp/. Let be a non-vanishing conormal to that cylinder at zf. After normalization, we get 

C“ = (- c9p ■ oo, tv) for some lo E S n 1 . Clearly, is spacelike. Let £ x » ! q p be the line on the cylinder 

C x e p through zo; then is conormal to it at zK 

To apply Proposition EU we claim that there is a lightlike line at z* normal to ^ so that this 
line is still in U. Suppose for a moment that this done. Then by by Proposition 13.11 we would get 
that / vanishes near that point, which would be a contradiction. Therefore, such a po would not 
exist, and in particular, / = 0 near £ Xl ,e 1 - 

To prove the claim, we are looking for a unit 9$ so that • (1 ,9$) = 0. This is equivalent to 
solving (9$ — c9) ■ cj = 0 for O'. It is easy to see that this is always possible to do and the solution 
closest to 9 is at its farthest distance from 9 when w = ± 9 ; then \9^ — c9\ = yj 1 — c 2 . This shows 
that 1 9^ — 9 1 < y/l — c 2 uniformly in aProperty (13.191) then proves the claim. □ 

4. The Lorentzian case 

4.1. Support theorems for analytic Lorentzian manifolds. Let ( M,g ) be a Lorentzian mani¬ 
fold now. Next theorem is an analog of Theorem l3.21 Since the global geometry of the null-geodesics 
in the Lorentzian case is non-trivial and in particular, one can have conjugate points, the assump¬ 
tions are stronger. 

Definition 4.1. Let S be a smooth surface near a point z E S and let F be a defining function 
so that S = f ?_ 1 (0) near z, d F(z) / 0, and declare {i 7 < 0}, to be the “interior” of M near z. 
Similarly, {F > 0} is the “exterior” of M near z. We say that S is strictly convex at z in the 
direction v € T Z S, ifX7 2 F(z)(v,v) > 0. 

We call S strictly lightlike-convex if it timelike, it is strictly convex at all lightlike (z,v) E TS, 
and every maximal lightlike geodesic tangent to S at some point has no other common points with 
S. 
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Here V 2 F is the Hessian of F, with V being the covariant derivative. This notion of convexity is 
equivalent to o 7 (s) < 0 for the geodesic 7 through x in the direction v\ and it is independent 
of the choice of F. 

Theorem 4.1. Let (. M,g ) be an analytic Lorentzian manifold. Let S be a timelike surface near 
a fixed point zq € S. Let 70 be a lightlike geodesic through zo tangent to S at zq. Assume that 
S is strictly convex at zq in the direction of 70 , and that k is analytic and non-vanishing near 
(>o,7oL)- Let f be a distribution, and let 70 be non-trapping in supp f. Let L K f{ 7 ) = 0 for all 
lightlike geodesics 7 near 70 . If f = 0 in the exterior of S near zq , then f = 0 near zo ■ 

Proof. By the Sato-Kawai-Kashiwara Theorem, see m and [28j, that we already used in the proof 
of Proposition 13.11 it is enough to prove that / is microlocally analytic in the direction of the 
conormal Co to S at zq. 

We follow the construction in Proposition 13.11 We can consider the former proof as a linearized 
version of the present one, when we replace S with its tangent plane at zq normal to Co> and the 
geodesics through zq by tangent lines. 

By the strict convexity assumption, for any lightlike geodesic 7 which is an 0(e), e <C 1, pertur¬ 
bation of 70 (in a fixed parameterization), the intersection of 7 with the interior of S, in any local 
chart has Euclidean length 0(y/e). This allows us to work in a fixed coordinate system (t, x) near zq. 
We choose a spacelike surface Sq through zo. We then choose semigeodesic coordinates (t,x) near 
zo = 0 normal to So, i.e., the lines (t, x) = (t, const.), are future pointing (the future direction being 
determined by 70 ) timelike geodesics normal to Sq and g is given locally by — df 2 + g a gdx a dx^, 
see [23]. Such coordinates are constructed by taking a normal field v to S normalized so that 
g{v,v) = — 1 and using it as initial directions of the geodesics x = const. We can arrange that 
zo = 0 and 70 ( 0 ) = ( 0 ,e„). 

In those coordinates, future pointing geodesics near 70 , close to So, are parameterized by their 
initial points x € S and the projection 6 of their tangents to TSq, i.e., Xx,e{s) is defined as the 
geodesic issued from ( 0 ,x) with 7 ^ 0 ( 0 ) = ( 1 ,6). 

For (x,9) close to (0, e n ), we then write 

L K f(x,0) = J n('y Xt B(s),j X} e{s))f('y X} e{s))ds. 

We chose 9 = 9(q) as in (|3.5p with |g| < e, where e is the number controlling the size of supp xn, 
see (13.101) . Then for 0 < e <C 1, A > 0, 

° = / e lXx< (XNL K f)(x,6(q))dx 

= Jf e lXiC< XN(x)K('y x ^ q) (s), 7 ^ 0 (g )(s))/( 7 a . i 9 ( q ) (s))ds dx. 

For every fixed q near q = 0, which fixes 9 = 9(q), the map (s,x) -A z = 7 ^ 0 ( 5 ) is a local 
diffeomorphism near zo = 0 by the Implicit Function Theorem. Let s\z,9) , x^(z, 9) be the inverse 
map. Since 7 x ,e(s) = (s,x + s9) + 0(s 2 ), we get the Taylor expansion 

(4.1) s\z,9(q)) = t + Oft 2 ), x\z, 9{q)) = x — t9(q) + 0(t 2 ), 

where z := (t,x). Those expansions can be justified by the Implicit Function Theorem. Make the 
change of variables (s, x) -A z above to get 

0 = J e lA0 XA i{x\z,9{q)))nJ{q,z)f{z) dz 


(4.2) 
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with (j)(z,£,q) = x$(z,6(q)) ■ £. Here, k is the weight in the new variables, and J is the related 
Jacobian. If g is Minkowski, we get <fi = (z — tO(q)) ■ £, which is the same function as in (13.121) . 

Set C = (<?, £)• Th en q = ( 0 , £ = (' = (Cl, ■ ■ ■ , Cn) and 

(4.3) <f>(z,C) = xHz,0(Co)) ■ ('■ 

Lemma 4.1. det <^(0, e n_1 ) = — 1. 

Proof. To compute (j> x ^(0,e n ~ 1 ), write first (recall that zo = t) 

(4.4) <falc=e"-i lZO =o = x i (z,e n )\ Z0=0 = z k , k = l,...,n. 

Therefore, 

^(O,^- 1 ) = 8i, k = 1,. ..,n, i = 0,1,... ,n. 

Therefore, det 4> z ({0, e n_1 ) = ^ Q 2 o(0, e n_1 ). One the other hand, the latter equals —1 as follows 
from (|4.1I) and (13.51) . □ 

We now get from (14.21) : 

(4.5) 0 = J jW z >Qa N (z, ()f(z) dz = 0 near C = e™" 1 , 

compare with (|3.12l) . with cln elliptic and analytic near (0, e"^ 1 ) but not for all (z, (). On the other 
hand, it satisfies a pseudo-analytic estimate of the type (13.101) . 

Similarly to the proof of Lemma 13.11 for w and q as (13.131) , in multiply (14.51) by the factor 

g iA(i(C—»?) 2 /2-0 Oj,C)) 

and integrate w.r.t. ( over the ball |£ — 77 1 < 5 with 0 < 6 <C e to get (13.141) with 

$ = C) - <l>(w, 0 + i(C - 7 ?) 2 / 2 - 

The rest of the proof follows closely those in [71133]. By Lemma [4. 11 |<3?g| has a lower bound outside 
any neighborhood of z = w for w localized as above and z in the support of the integrand. This 
allows us to integrate by parts to get ()3.151) in this case and choose N ~ A /(Ce) to make the r.h.s. 
of ([3.151) exponentially small with A. The phase function has an analytic extension for ( in some 
complex neighborhood of £0 = e n-1 . By Lemma 14.11 <f^{zX) = ^(to, C) f° r such £ and z, w close 
to 0 implies z = w. Therefore, the critical point ( c = q of <h w.r.t. ( is real only when z = w and 
at that point, > 0 ; and it is unique and complex otherwise, still satisfying that inequality by 

a perturbation argument, when 0 < <5 <S e. Then we get (13.171) with if defined as in (13.161) . Then 
we conclude as in mm that (0,6) £ WF a (/). 

This arguments so far work if / is a continuous function, for example. If / is a distribution, 
as stated, we need to take a smooth cutoff xs and consider the z-integrals above in distribution 
sense. □ 

Proof of Theorem \2.2l Assume that the statement of the theorem is not true. Let do € [0,1) be 
the infimum of all a for which i ?_ 1 (cr) D supp / 7 ^ 0. Then / = 0 in the “exterior” i ?_ 1 (0, 00 ) of 
Sq = F^ 1 (ao), and supp / has a common zq point with Sq. The latter follows from a compactness 
argument. In particular, co > 0. Since So is timelike by assumption, there is a lightlike geodesic 70 
through zo tangent to So which does not hit So again by the strict convexity assumption. Near zo, 
the geodesic 70 lies in the exterior I ?_ 1 [ 0 , < 7 o] by the local part of the strict convexity assumption; 
and this is also true globally by the global part of that assumption. By Theorem 14.11 / = 0 near 
every common point of So and supp/, which is a contradiction. □ 
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